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ABSTRACT 



Context. Full Stokes filter-polarimeters are key instruments for investigating the rapid evolution of magnetic structures 
on the solar surface. To this end, the image quality is routinely improved using a-posteriori image reconstruction 
methods. 

Aims. We analyze the robustness of circular polarization asymmetries to phase- diversity image reconstruction techniques. 
Methods. We use snapshots of magneto-hydrodynamical simulations carried out with different initial conditions to 
synthesize spectra of the magnetically sensitive Fe i line at 5250.2 A. We degrade the synthetic profiles spatially and 
spectrally to simulate observations with the IMaX full Stokes filter-polarimeter. We also simulate the focused/defocused 
pairs of images used by the phase- diversity algorithm for reconstruction and the polarimetric modulation scheme. We 
assume that standard optimization methods are able to infer the projection of the wavefront on the Zernike polynomials 
with 10% precision. We also consider the less favorable case of 25% precision. We obtain reconstructed monochromatic 
modulated images that are later demodulated and compared with the original maps. 

Results. Although asymmetries are often difficult to define in the quiet Sun due to the complexity of the Stokes V 
profiles, we show how asymmetries are degraded with spatial and spectral smearing. The results indicate that, although 
image reconstruction techniques reduce the spatial smearing, they can modify the asymmetries of the profiles, mainly 
caused by the appearance of spatially-correlated noise. 
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1. Introduction 

The observation of Stokes profiles in spectral lines are a 
valuable tool to infer information about the thermodynam- 
ical and magnetic properties of the solar plasma. This in- 
formation is encoded in the amplitude and the shape of the 
Stokes profiles. Therefore, it is important to avoid any ef- 
fect that perturbs the shape because it can crucially modify 
the information encoded in the profile. 

Of special relevance is the asymmetry of the circular po- 
larization profile, which is known to be related to the corre- 
lation between velocity and magnetic field gradients along 
the line-of-sight (LOS; Illing et all 119751 ). This effect has 
been exploited to build atmospheric models with such gra- 
dients to explain asymmetries around magnetic flux concen- 
tra tions fSolanki & Pahlke "l988"; Grossmann-Doerth et al. 
[1988: Sanchez Almeida et al. 1989; Solanki & Montavon 
Il99 3). The situation is especially relevant in the weakly- 
magnetized zones of the quiet Sun away from active re- 
gions, where Stokes V profiles prese nt a variety of shapes 
with strongly a symmetric profiles (Sigwart h et al.l 
Sanchez Almeida fc Lites 2000; Khomenk o et al.l 



Socas-Navarro et an 2004t IM artmez Gonzal ez et al.l 



1999; 



2003; 



20081; 



Viticchie et al. I l201l[lviticcliie fc Sanchez Almeidall201lD . 



Earth-based observations are always affected by the 
disturbing effect of the atmosphere. As a consequence, 
the diffraction limit of the telescope is practically never 
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reached. It is often impossible to overcome the 1" limit 
if the observations are not accompanied by an adaptive 
optics syste m and powerful p o st-proc essing methods. For 
this reason, iKhomenko et al.l (|2005b[ ) and IShelvag et al.l 
(2007) analyzed the effect of spatial smearing on the Stokes 
asymmetries observed in the pair of Fe I lines at 630 nm 
and at 1565 nm. They concluded that asymmetries are 
heavily disturbed by the lack of spatial resolution. They 
even discovered that it is possible to find regions in which 
the field polarity is different i n the two s pectral regions 
(Sa nchez Almeida et aD 120031 : iRezaei et al.l 120071 ). some- 
thing fully attributed to the lack of spatial resolution. The 
fundamental reason is that the shape of Stokes V profiles 
in weakly-magnetized regions are very complex and they 
change in scales much smaller than the resolution element 
of the largest Earth-based telescopes even in the hypothet- 
ical absence of atmosphere. 

Even if t he telescope is put in a balloon at 40 km height 
like Sunrise ([Solanki et al.ll201Ql ), there is some remaining 
atmosphere (albeit small) which, together with the intrin- 
sic aberrations of the telescope and instruments, can mod- 
ify the observations. For this reason, the IMaX instrument 
onboard Sunrise was designed to us e the phase-diversity 
post-facto reconstruction algorithm (|Paxman et al.l 119921 
'Vargas Dommguez"2009'). In this paper, we focus on the 
interesting problem of testing how the spatial and spectral 
degradation produced by the Sunrise/IMaX combination 
affects circular polarization asymmetries and if the post- 
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reconstruction algorithms can help us extract reliable in- 
formation (comparable to the unperturbed synthetic case) 
about them from degraded data. To this end, we use pro- 
files synthesized on 3D models of the solar photosphere. We 
degrade them to simulate the observational conditions of 
Sunrise/IMaX and reconstruct the images using the phase- 
diversity algorithm. 



2. Data processing 

2.1. Simulations and spectral synthesis 

The snapshots that we have used correspond to individ- 
ual time steps of a 3D magneto-hydrodynamical simulation 
of solar magn eto-convection don e with the MURAM code 
(IVQglerll2QQ3l : IVoder et al.ll2QQ5l : ICameron et al.ll2Qlll) . An 
initially vertical magnetic field of 200 G strength was intro- 
duced into already developed purely hydrodynamical con- 
vection. The simulation box was split into 4 parts with 
the opposite polarities of the magnetic field in the adjacent 
parts. This magnetic field evolved self-consistently with 
convective motions. The redistribution of the magnetic field 
led to almost exponential decrease with time of its average 
unsigned value over the simulation domain. The snapshots 
used in the present work were taken 17, 36 and 112 min- 
utes after the magnetic field was introduced. At these time 
moments, the average unsigned magnetic field strength at 
photospheric base was of 140 G, 80 G and 30 G. The last 
snapshot is representative of quiet Sun internetwork be- 
cause i t is close to the value obtained bv iKhomenko et al.l 
(|2005at ) , while the first two are representative of bipolar en- 
hanced network. The computational box has 288 x 288 x 100 
grid points, with the largest number of points in the hori- 
zontal directions. The grid step in the horizontal direction 
is roughly 20 km, with a total horizontal extent of 6 x 6 
Mm^. The grid size in the vertical direction is ~14 km, 
with an ensuing total height of 1400 km. On average, the 
spatial extent above and below optical depth unity at 5000 
A amounts to 600 km and 800 km, respectively. 

We have computed the Stokes spectra of the Fe 
I 5250.2089 A line formed at solar disc centre (/i = 
1). The calculation was done u sing the SIR code 
(iRmz Cobo fc del Toro Inii^ fTool ) in LTE approxima- 
tion. The spectra were computed in a 2.5D approach, i.e. 
individually for each vertical column of the snapshots, cor- 
responding to the LOS direction. We used the following 
atomic parameters: lower level energy Ei = 0.121 eV, 
\og gf = —4.96 and iron abundance Afq=7.50. The detailed 
description of this particular simulation run and analysis of 
Stokes spe ctra of th is and other Fe I lines can be found 
in Khomen ko et al.l (j2b05b) and iKhomenko fc ColladosI 
(1200711 . 

Since we simulate the effect of image degradation on 
the polarimetric measurements, we first have to modulate 
the resulting monochromatic synthetic images of the Stokes 
parameters. These modulated images will be distorted by 
the presence of aberrations and measured on the CCD. For 
simplicity, we apply the ideal modulation scheme: 

1111" 
1 1-1-1 
1-1-1 1 
1-1 1-1 
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1 
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Distance [Mm] 



and neglect any cross- talk between Stokes parameters. 
Differences between the real modulation matrix for IMaX 



Fig. 1. Continuum (left panels) and total polarization 
(right panels) computed as / (Q^ + t/^ + V^)dA/ / IdX for 
the three snapshots considered of 30 G (upper panels), 80 G 
(middle panels) and 140 G (lower panels) of average mag- 
netic field at logrsooo = — 1- 



and the ideal one presented above will fundamentally af- 
fect the efficiency of the modulation scheme. Consequently, 
the only difference lies on the final noise level. The re- 
sulting modulated images are spatially degraded, noise is 
added and the reconstructed image is obtained through 
standard phase-diversity formulae. The ensuing Stokes pa- 
rameters are obtained by applying the demodulation matrix 
D = O-^. 



2.2. Degradation and phase-diversity reconstruction 

The synthetic profiles at full resolution are degraded to 
the IMaX spatial and spectral resolution. IMaX is a filter- 
polarimetric instrument based on Fabry-Perot interferom- 
eters, so that the spectral point spread function (PSF) of 
the double passage is known analytically. This PSF has very 
extended tails that would make it necessary to synthesize a 
large portion of the spectrum to correctly take i nto account 
the sm earing. Instead, we follow the approach of lLagg et al.l 
(j2010[ ) who substituted the real PSF by a Gaussian PSF. 
It has been verified empirically that a Gaussian of 85 mA 
(- 2.9 km s"^) of full width at half-maximum (FWHM) 
gives a good representation of the effect of the real PSF, 
taking into account the contribution of secondary lobes. 
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The smeared profiles are then sampled at 12 wavelength 
points to simulate the LI 2-2 observing mode of IMaX 
(|Martmez Pillet et alJ[2QTl h. from -192.5 mA to +192.5 
mA in steps of 35 mA. 

Concerning the spatial smearing, we have simulated all 
image degradation effects that have been considered for 
Sunrise/IMaX f see ' Vargas D ominguezl l2QQ9l for a detailed 
analysis). Adding all aberrations, a value of A/5 for the 
rms wavefront error (WFE) has bee n measured as the ac- 
tual performance of the instrument (|Martfnez Pillet et al.l 
l2Qlll ) . We consider the same aberration but divided in sev- 
eral contributions. The considered telescopic aberrations 
(up to the 45th Zernike polynomial) produce a wavefront 
with an rms WFE of A/9. These aberrations are slowly 
varying in time so we consider them fixed during the full po- 
larimetric modulation and spectral scanning. The value of 
the coefficients associated to each Zernike polynomial that 
conform the wavefront telescopic aberration are described 
by uncorrelated Gaussian distributions. The atmosphere 
is accounted for by considering a wavefron t with turbu- 
lent K olmogorov statistics (Noll 1976; Lofdah Ffc Scharmerl 
Il994t ). again with an rms WFE of A/9. While this atmo- 
spheric contribution might be considered large for the case 
of IMaX/ SUNRISE, including it here allows our results to 
be valid for ground-based observations with adaptive op- 
tics systems. This atmospheric contribution might be even 
considered large for the case of IMaX/SUNRISE. Assuming 
the presence of an image stabilization system, we neglect 
the presence of tip and tilt. We also consider a ripple WFE 
of A/9 associated to the polishing that we simulate as a 
screen with a von Karman spectrum with an external scale 
of 30 cm that simulates the diameter of the polishing tool. 
Considering the previous contributions, the total rms of the 
WFE obtained as the addition in quadrature is of the order 
of A/5. Furthermore, we also take into account the effect of 
the finite size of the CCD pixels, which do not sample the 
image with Dirac delta functions but with a small amount 
of diffraction. Finally, the resulting PSF is then convolved 
with the monochromatic modulated images and noise at the 
level of 10~^ in units of the continuum intensity is added. 
After demodulation, we obtain the degraded monochro- 
matic images of the Stokes parameters. Additionally, since 
we apply image reconstruction techniques based on the 
phase-diversity technique, we also simulate images with ex- 
actly the same wavefront but including an additional opti- 
mal defocus, in which the peak-to- yalley defocus w avefront 
equals the observation wavelength (|Lee et al.lll999f ) . For the 
case of IMaX, this optimal defocus amounts to ~ 1.81A. 

Our aim is to compare the asymmetries of Stokes V 
in the degraded profiles with those found in the original 
one and in the phase-diversity reconstructed data. Since 
we know exactly the PSFs of the focused and defocused im- 
ages, we could apply a standard phase-diversity recon struc- 
tion algorithms (Gonsalves fc Chidlaw 1979 ; Paxman et al.l 
1 19921 : iLofdahl fc Scharrnerl Il994l ) and obtain a perfect re- 
construction using: 



(2) 



where Pq, Pk^ Dq and Dk are the Fourier transforms of 
the focused PSF, defocused PSF, focused image and defo- 
cused image, respectively, and F is the Fourier transform 
of the reconstructed image. In real situations, the PSFs are 



not known and the reconstruction algorithm is preceded 
by an iterative process in which the PSFs are estimated 
by maximizing a merit function. In order to overcome this 
iterative process but still consider uncertainties in the esti- 
mation of the PSFs, we make the assumption that the pro- 
jection of the final wavefront on every Zernike polynomial is 
known with 10% relative error. This figure is representative 
of extensive simulations carried out with IMaX data (e.g., 
IVargas Dommguezll2QQ9[ ). This error is representing the un- 
certainties that one would have when using optimization 
methods for the estimation of the wavefront. From the re- 
constructed focused and defocused PS Fs, we obtain the re - 
constructed ima^e fol lowing Gonsalves fc Chidlawl (119791): 
iPaxman et HI (|1992[ ): ILofdahl fc Scharmerl (11991 )7 Because 



of the presence of noise and the diffraction limit of the tele- 
scope, it is fundament al to filter out high freq uencies from 
the resulting image (Lofd ahl fc ScharmedfTggl . These high 
frequencies cannot be recovered. Additionally, the noise on 
the reconstructed images present a Gaussian distribution 
with a variance larger than that of the original focused and 
defocused images. Th e main probleni is th at this noise is 
spatially correlated (Mevnadi er et al.lll999f ). This leads to 
an aesthetically unappealing web-like noise which avoids us- 
ing many profiles in the quiet- Sun because their amplitudes 
are already very low. Furthermore, this noise introduces ar- 
tifacts that perturb the appear ance of the Stoke s p rofiles . 
This effect is also discussed by IPuschmann fc Bec3 (|2Qllf ) 
with different deconvolution methods. 



3. Circular polarization asymmetries 

Circular polarization asymmetries are usually defined in 
terms of the relative area {SA) and amplitude (Sa) imbal- 
ances bet ween the blue and the red lobes of Stokes V pro- 
files fe.g. JSolanki fc Stenfiol[l986 ): 



Sa 
SA 



Clh — Clr 



= S 



JV{X)dX 
I\V{X)\dX' 



(3) 
(4) 



where and refer to the amplitude of the red and blue 
lobes, respectively. The factor s gives the sign of the area 
asymmetry, and it is chosen equal to the sign of the blu est 
peak of Stokes V, following (i Martmez Pillet et al.|[T997l ). 

These definitions of asymmetries can only be strictly 
applied to standard two-lobed antisymmetric profiles (one- 
lobed profiles can also be considered as a degenerate two- 
lobed profile). However, simulated Stokes V profiles of- 
ten show a panoply of profiles whose spectral shapes are 
completely different to the standard antisymmetric profile. 
Table [T] displays the percentage of Stokes V profiles with 
different number of lobes on the three snapshots consid- 
ered. Although the majority of profiles are standard two- 
lobed profiles, there are many locations where 1-lobed or 
3-lobed profiles are found. This situation is more critical 
for the case in which the full spatial and spectral resolution 
is considered. To overcome this problem, we have restricted 
our analysis to one-lobed and two-lobed profiles, leaving for 
the future the more complex analysis of profiles with more 
lobes. In such a case, we think that it will be necessary to 
update the definition of asymmetries. Interestingly, recent 
high spatial resolution observations (Martinez Gonzalez et 
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Table 1. Statistical information about profiles 



Snapshot 


Case 


Above 


1 lobe 


2 lobes 


3 lobes 


4 lobes 


Amplitude 


Area 






threshold"" 










assymetry^ 


asymmetry^ 


30 G 


Fuh (IMaX) 


36% 


14.61% 


80.36% 


4.32% 


0.05% 


8.65 =b 38.6 


6.68 ± 42.4 


- 


Degraded (IMaX) 


11% 


0.75% 


99.25% 


0.01% 


0.01% 


8.65 zb 38.6 


-5.72 ± 15.8 


- 


Reconstructed (IMaX) 


5% 


17.44% 


82.48% 


0.10% 


0.03% 


1.79 =b 14.7 


15.58 ± 34.92 


- 


Full (all A) 


54% 


25.12% 


64.78% 


9.88% 


0.07% 


28.10 ±33.9 


3.02 ± 29.2 


80 G 


Fuh (IMaX) 


68% 


14.00% 


80.90% 


4.90% 


0.04% 


9.04 ± 39.2 


8.12 ± 42.5 




Degraded (IMaX) 


34% 


1.27% 


98.71% 


0.02% 


0.01% 


4.53 =b 16.2 


-0.69 ± 18.0 




Reconstructed (IMaX) 


21% 


14.60% 


84.74% 


0.65% 


0.01% 


22.10 ±31.8 


16.54 ±32.6 




Full (ah A) 


85% 


16.26% 


69.55% 


14.07% 


0.06% 


27.67 ± 37.6 


5.00 ± 30.6 


140 G 


Fuh (IMaX) 


85% 


12.61% 


84.30% 


3.05% 


0.01% 


18.07 ±40.0 


19.12 ±43.5 




Degraded (IMaX) 


71% 


2.42% 


97.32% 


0.26% 


0.00% 


13.53 ± 20.6 


11.46 ±23.4 




Reconstructed (IMaX) 


34% 


13.19 


85.49% 


1.31% 


0.00% 


21.66 ±32.3 


20.18 ±34.2 




Full (all A) 


95% 


7.04% 


80.55% 


12.34% 


0.05% 


28.87 ± 38.4 


13.50 ±32.2 



^ Percentage of profiles with Stokes V signal above the noise threshold. 
^ Mean and standard deviation of amplitude asymmetry. 
^ Mean and standard deviation of area asymmetry. 



al. 2011, in preparation) also show the appearance of an 
increasingly larger fraction of non-standard profiles. 

In order to analyze asymmetries, we consider four dif- 
ferent cases that will help us understand how Stokes V 
profiles are distorted by the presence of spatial and spec- 
tral smearing. The first case, which we label as "Full (all 
A)" consists of the line profiles at full spatial and spec- 
tral resolution. The next, labeled as "Full (IMaX)", are 
the profiles at full spatial resolution but degraded to the 
IMaX LI 2-2 observing mode spectral resolution. We also 
consider two observationally realistic cases. The first one 
is the fully degraded, labeled as "Degraded (IMaX)". The 
standard deviation of the noise added a-posteriori in units 
of the continuum intensity in these snapshots is 10~^. The 
last one is the phase-diversity reconstructed case, labeled 
"Reconstructed (IMaX)". The standard deviation of the 
noise increases to 3 — 4 x 10~^ in units of the continuum 
intensity after the reconstruction. 

3.1. Spatial variation 

Figures [2] and [3] show the amplitude and area asymmetries 
for all considered cases, respectively. The upper panels show 
all map locations where the Stokes V profiles are above 
the threshold for each individual case. For the sake of a 
better visualization, the bottom panels display only map 
locations whose Stokes V profiles are above 3 times the 
value of the standard deviation of the noise in all cases 
simultaneously. In this case, the "Reconstructed (IMaX)" 
case is the one controhing such map locations because of 
the enhanced noise level. 

From left to right we show what happens when the spa- 
tial/spectral resolution is degraded, while the last panel 
displays the phase-diversity reconstructed case. The three 
cases considered (each row) correspond to the three snap- 
shots with different seed magnetic fields, where only profiles 
above the threshold of 3 times the noise level have been con- 
sidered. Obviously, the number of selected profiles increases 
as the seed field increases because of the enhanced Stokes 
V amplitudes. 

Concerning the high-resolution maps, we find predom- 
inantly positive area asy mmetries in g^ranu le s, simi lar to 
the previous analysis of iKhomenko et al.l (|2005b[ ) and 



IShelvag et aD (|2007f ) for the Fe I lines at 630 nm and 1.56 
/j,m. This result is compatible with the predominance of 
posit ive area asymmetries in gra nules found in observations 
(e.g. JSocas-N avarro et al.]|2004 ^. This predominance is also 
found for amplitude asymmetries. Negative asymmetries 
are found in intergranular lanes and in strong kG acum- 
mulations. This result found in the high-resolution maps 
is essentially incompat i ble w ith the observational results 
of 'Socas-Nav arro et al.l (|2004l ), who find positive area and 
amplitude asymmetries in lanes and the network. We show 
below that both results can be reconciled after image degra- 
dation. 

Asymmetries in the high-resolution maps can be under- 
stood under the fiux tube model with a canopy developing 
with height (|Bellot Rubio et al.l l2000[). The inner parts of 
the flux tubes coincide with intergranular lanes in the sim- 
ulations and present strongly negative area asymmetries 
with not so strongly negative amplitude asymmetries. They 
are generated by strong downflows a ssociated with a fiel d 
strength that decreases with height ([Shelvag et al.l [20071 ). 
The granule-intergranule transition with positive area and 
amplitude asymmetries can be explained by rays crossing 
the canopy of the flux tube. 

Apparently, the spectral degradation of the data at full 
spatial resolution has a relatively small impact on the spa- 
tial distribution of asymmetries, as seen from a comparison 
of the first and second columns. The filamentary structure 
of the asymmetries associated with intergranular lanes is 
maintained, while the profiles observed in granular regions 
are lost because the spectral degradation pushes them be- 
low the noise threshold. However, we can see that the spec- 
trally degraded data presents larger (in absolute value) area 
and amplitude asymmetries than the original dataset, being 
this difference especially relevant on the cores of intergran- 
ular lanes. This can be easily understood because spectral 
degradation tends t o increase relative area and amplitude 
asvmmetries Solanki fc Stenfiol ((1986), although the depen- 
dency on the specific shape of the Stokes V profile is crucial. 

Spatial degradation changes completely the scenario. 
The main reason is that spatial smearing produces two ef- 
fects: i) profiles with very different asymmetries are mixed, 
and ii) a decrease on the Stokes V amplitude, thus push- 
ing the signal in many pixels below the noise threshold. 
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Fig. 2. Amplitude asymmetry maps for the four considered smearing cases (columns) and for the three values of the seed 
magnetic field used (rows): 30 G, 80 G and 140 G. The upper panels show the map locations where the Stokes V profiles 
are above the threshold in each individual case. The lower panels display only locations where the Stokes V profiles are 
above the threshold simultaneously in the four considered cases. The small green lines at positions (2,1) for the upper 
panels, (0.5,0.5) for the middle panels and (2.3,3.8) for the lower panels show the location of the cuts along which the 
Stokes V profiles of Fig. |4] have been obtained. 
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Still, the pixels associated with the strongest signals re- 
main (11%, 38% and 70% of the pixels for the 30 G, 80 
G and 140 G snapshots, respectively) and the filamentary 
structure is somehow lost and transformed into a blobby 



structure. Typically, asymmetries tend to be organized in- 
side these blobs so that zero or very small asymmetries are 
found in the cores while they increase in absolute value to- 
wards the borders, where profiles are lost below the noise 
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Fig. 4. Stokes V profiles normalized to the peak amplitude across the cuts shown in Fig. [2] and [3] in green. Black profiles 
correspond to the "Full (IMaX)" case, while green ones are associated with the "Full (all A)" case. Red profiles correspond 
to the "Degraded (IMaX)" case and blue profiles are the "Reconstructed (IMaX)" ones. 



threshold. This behavior is not exactly the same of what 
happens at larger spatial resolutions, where the cores of in- 
tergranular lanes present negative asymmetries while they 
become positive on the granule-intergranule transition re- 
gion. From a general perspective, we can say that spatial 
smearing makes the negative asymmetry of narrow inter- 
granular lanes appear with smaller absolute values while 
broadening the region of positive asymmetry. 

The area and amplitude asymmetries found in the spa- 
tially degraded snapshots tend to be more compatible with 
observations. On average, asymmetries tend to be positive, 
with a larger absolute value for amplitude asymmetries 
than for area asymmetries. Stokes V profiles in integranular 
lanes have now negative area asymmetries not larger than 
15-20% or even compatible with zero. The only exceptions 
are isolated patches with strong negative asymmetries. 

It is evident from Figs. [2] and [3] that area asymmetries 
are slightly less affected by the spatial and spectral smear- 
ing than amplitude asymmetries. It is specially relevant in 
the central parts of internetwork lanes, where strong neg- 
ative amplitude asymmetries show up when reducing the 
spectral resolution and these are transformed into almost 
amplitude symmetric profiles when degrading the spatial 
resolution. On the contrary, relevant negative area asymme- 



tries are still present on internetwork lanes in the spatially 
and spectrally smeared data. 

The main point of this work is to verify to what extend 
asymmetries are recovered after image reconstruction tech- 
niques. The rightmost panels of Figs. [2] and [3] demonstrate 
that an improvement on the spatial location of asymme- 
tries occurs. Obviously, the enhanced noise level and its 
spatially correlated character produces that the number of 
pixels with detected signals decreases (6%, 22% and 43% of 
the pixels for the 30 G, 80 G and 140 G snapshots, respec- 
tively). Comparing with Fig. [TJ only pixels associated to 
strong magnetic concentrations can be correctly analyzed. 
Consequently, one should be careful when analyzing recon- 
structed data to take into account this bias. Several conclu- 
sions can be extracted from the reconstructed images. First, 
the structures are much more compact than in the degraded 
case, a consequence of the efficient reduction of aberrations 
performed by the phase-diversity algorithm, even with 10% 
relative errors on the projections of the wavefront on the 
Zernike polynomials. Second, the large gradient of asymme- 
tries from the core of intergranular lanes to the surround- 
ings is partially restored. However, the large negative value 
of SA and 6a found in the very central cores are not fully 
recovered. Generally, in these regions, image reconstruction 
cannot generate strong negative asymmetries if they were 
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Fig. 5. Left panels: area asymmetries obtained on the degraded snapshots (black for "Full (IMaX)" case, red for 
"Degraded (IMaX)" and blue for "Reconstructed (IMaX)") versus the original "Full (all A)" area asymmetries. Right 
panels: same but for amplitude asymmetry. The upper panels show results when the wavefront is known with 10% 
uncertainty, while the lower panels display what happens when the uncertainty increases to 25%. 



not present in the degraded images. Likewise, if they were 
present in the degraded maps, they will be enhanced in the 
image reconstruction process. 

3.2. Line profiles 

Synthetic Stokes V profiles for the three snapshots consid- 
ered are shown in Fig. |4] for the cuts indicated in Figs. [2] 
and [3] in green. The Stokes profiles are normalized to peak 
amplitude so one has to take into account that spatially 
and/or spectrally averaged signals possess a smaller ampli- 
tude. The profiles of the full spectral and spatial resolution 
are shown in green. The profiles of the full spatial resolu- 
tion but spectral resolution degraded to IMaX are shown 
in black. The fully degraded profiles are shown in red and 
the phase-diversity recovered ones are plotted in blue. 

The first impression is that the full spatial and spec- 
tral resolution profiles are very narrow and the presence of 
several components (several lobes on the red lobe) along 
the LOS can be easily witnessed. When the spectral res- 
olution is degraded to 85 mA, the profiles are broadened, 
with strong modifications of the relative amplitudes of the 
two lobes. As a byproduct, the full spatially /spectrally de- 
graded profiles tend to be much more symmetric and only in 
those cases in which the full resolution profile is extremely 



asymmetric, some asymmetry remains on the degraded pro- 
files. 

It is interesting to note that the phase-diversity recon- 
structed profiles tend towards the full spatial resolution case 
with IMaX spectral resolution, as expected. This means 
that an image reconstruction algorithm is able to partially 
extract the information necessary to recover the Stokes V 
profiles. This is especially relevant for the very asymmetric 
profiles. 

3.3. Statistical properties 

We have tested that image reconstruction based on phase- 
diversity does a good job at cancelling the spatial smearing 
of polar imetric signals. Furthermore, although many points 
are lost below the enhanced noise level, the spatial appear- 
ance of asymmetries is not strongly modified. However, it 
is important to verify other properties. We focus now on 
the statistical properties of asymmetries. 

Fig. [5] displays the area (left panels) and amplitude 
(right panels) asymmetries computed on the degraded 
snapshots versus those obtained on the full resolution snap- 
shots. The upper panels correspond to the case in which 
the wavefront is known with 10% uncertainty, whereas the 
lower panels are associated to the case in which the un- 
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certainty increases to 25%. Red points correspond to the 
"Degraded (IMaX)" case. Black points correspond to the 
"Full (IMaX)" case, while the blue points are associated 
to the "Reconstructed (IMaX)" case. Since the number 
of points with profiles above the noise level is different 
for each case, we consider only those points which pro- 
duce signal above the threshold in the four cases simul- 
taneously. Essentially, this leaves only those points on the 
phase-diversity reconstructed images, as shown in the lower 
panels of Figs. [2] and [3l 

The spectral degradation produced by IMaX induces 
that area asymmetries (in absolute va lue) a re systemat- 
ically overestimated Solan ki fc St e nflol (|l986h because the 
black points are above (for positive asymmetries) and below 
(for negative asymmetries) the main diagonal. Concerning 
the amplitude asymmetry at IMaX spectral resolution, the 
lack of spectral resolution induces a spread around the 
correct value. When the spatial resolution is degraded to 
IMaX resolution, we find that both area and amplitude 
asymmetries are always underestimated with respect to the 
original one. However, while positive asymmetries are un- 
derestimated, negative asymmetries are severely lost. The 
reason is that these points with negative asymmetries are 
associated to profiles with low Stokes V amplitude and 
are quickly lost below the noise level (see iKhomenko et al.l 
l2QQ5b: Shelyag et al. 2007, for examples on simulations). 

The phase-diversity reconstruction leads to a correction 
of the area and amplitude asymmetries such that the over- 
lap with the black points is enhanced. This overlapping is 
evident for the case in which the wavefront is known with 
10% accuracy but less evident when the uncertainty raises 
to 25%. As a byproduct, some portion of the large negative 
asymmetries appearing on the central parts of intergranular 
lanes (specifically, those surrounded by large positive asym- 
metries) are recovered. On the contrary, isolated strong neg- 
ative asymmetries (like that found at position (2.5,3) on the 
upper panel of Figs. [2] and [3]) are lost below the noise level. 
Consequently, we conclude that the image reconstruction 
algorithm is able to partially recover the main properties 
of area and amplitude asymmetries. 

Interestingly, area asymmetries can be corrected to be 
close to the original ones if they are divided by a factor 
between 0.3 and 0.5. This figure is slightly dependent on 
the mean magnetic field of the simulation. A correction 
factor of ~ 0.4 is a good compromise, although the spatial 
smearing is still present and strong asymmetries are fully 
lost. 

4. Conclusions 

We have analyzed the impact of phase-diversity image re- 
construction algorithms on the Stokes V asymmetries. To 
this end, we have employed three snapshots of 3D MHD 
models of solar magneto-convection with different mean 
magnetic fields, as representative of regions with different 
magnetization on the solar surface. We have synthesized 
the Fe I line at 5250.2089 A introduced spatial and spec- 
tral smearing to simulate observations with IMaX, and ap- 
plied a phase-diversity reconstruction assuming imperfect 
knowledge of the wavefront. We conclude that asymmetries 
are modified by the smearing process and that this effect 
can be alleviated using the reconstruction algorithm. Area 
and amplitude asymmetries tend to be overestimated af- 
ter the spectral smearing. This is compensated by the spa- 



tial smearing, which tends to underestimate all asymme- 
tries. Large negative asymmetries are lost below the noise 
level because they are usually associated to low-amplitude 
Stokes V profiles. Reconstruction successfully reduces the 
underestimation of asymmetries and partially recovers the 
main characteristics. However, patches of isolated large neg- 
ative asymmetries cannot be recovered. We point out that 
a large fraction of the pixels cannot be analyzed in the re- 
constructed data because of the enhanced noise. This noise 
has spatial correlation, which introduces additional prob- 
lems. The individual line profiles tend to look more similar 
to those of the full spatial resolution case with degraded 
spectral resolution. 
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